Rapid directional alignment of velocity and magnetic field in magnetohydrodynamic 

turbulence 
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We show that local directional alignment of the velocity and magnetic field fluctuations occurs 
rapidly in magnetohydrodynamics for a variety of parameters. This is observed both in direct 
numerical simulations and in solar wind data. The phenomenon is due to an alignment between 
the magnetic field and either pressure gradients or shear-associated kinetic energy gradients. A 
similar alignment, of velocity and vorticity, occurs in the Navier Stokes fluid case. This may be the 
most rapid and robust relaxation process in turbulent flows, and leads to a local weakening of the 
nonlinear terms in the small scale vorticity and current structures where alignment takes place. 



In magnetohydrodynamic (MHD) turbulence, the fluc- 
tuating magnetic field b and velocity fluctuation v enter 
on nearly equal footing. One consequence is that the non- 
linear MHD equations are in effect linearized when the 
fluctuation components of the magnetic field (in Alfvcn 
speed units) are everywhere equal (or opposite) to the 
velocity field. Such "Alfvenic" states have long been 
thought to be favored in relaxation processes Near 
Alfvenic states are observed in the solar wind plasma Q , 
mostly in the inner heliosphere Q ■ Global evolution to- 
wards the Alfvenic state, or "dynamic alignment," 
when it occurs, requires many nonlinear eddy turnover 
times. Here we describe a related, more rapid relaxation 
process, in which local, directional near-alignment of v 
and b emerges in less than one turnover time, for a wide 
variety of turbulence parameters. This process need not 
be associated with global alignment, but rather occurs 
independently in numerous localized patches. 

Dynamic alignment competes with other MHD relax- 
ation processes as shown in simulations @ where, for 
some parameters, it does not occur, or is incompletely re- 
alized. Solar wind observations also show that the degree 
of Alfvcnicity tends to decrease with increasing heliocen- 
tric distance, in apparent contradiction to the dynamic 
alignment principle. There arc suggestions that direc- 
tional alignment (a necessary condition for global dy- 
namic alignment) may be more ubiquitous. When MHD 
relaxation is described by a constrained energy mini- 
mization principle 0, cross helicity (Alfvenicity) and a 
magnetic invariant (helicity in three dimensions - 3D; or 
mean square flux function in two dimensions - 2D) are 
held constant, while energy is minimized. The emergent 
Euler-Lagrangc equations predict final states, in both 2D 
and 3D, and for all parameters, in which v and b are 
directionally aligned or anti-aligned. This theory is rea- 
sonably well confirmed by direct numerical simulations. 
Evidently, in the long time limit for decaying MHD tur- 
bulence, pointwise directional alignment is obtained more 
easily than is the global Alfvenic state. Below we show, 
using MHD numerical experiments, that local directional 



alignment is even more robust, occurs more rapidly, and 
appears locally, in patches. The distribution of alignment 
angle in the solar wind is shown to be consistent with this 
picture. This rapid relaxation can be understood by an 
elementary examination of the MHD equations. While 
these features appear not to have been fully recognized 
previously, the situation is in fact analogous to the local 
emergence of Beltrami flows Q in hydrodynamics. 

MHD and Alfvenic states. In familiar dimensionless 
( "Alfvenic" ) units, the equations of incompressible MHD 
are 

— +v-Vv = -VP+jxb + i/V 2 v (1) 
ot 



b ■ Vv + r/V 2 b, 



(2) 



with V • v = V • b = 0. Here, v is the velocity field, and 
b is the magnetic field, related to the electric current 
density jbyVxb=j;'Pis the pressure. The viscosity 
v and magnetic diffusivity i] define mechanical and mag- 
netic Reynolds numbers respectively as Ry = LU/u and 

R M = LU/r). Here U = (|v| 2 ) 1/2 , with (...) denoting a 
spatial average, and L is a length scale associated with 
the large-scale flow (integral length scale) . The total en- 
ergy E = E v + E b = (|v| 2 + |b| 2 )/2, the cross helicity 
H c = (v ■ b), and the magnetic helicity H m = (a • b) 
are ideal (v = rj = 0) invariants in 3D. Here b = V x a. 
Dimensionless measures of the helicities are a c = 2H C /E 
and a m = (El — Eh) / 'Eb, where El and En are magnetic 
energy in left- and right-handed magnetic polarizations, 
respectively. 

Simulations. We consider several sets of simulations 
(see Table I), in which the MHD equations are integrated 
in a spatially periodic domain of side 2tt, using a second 
order Runge-Kutta method, and either 2/3-rule dealiased 
H, or non-dealiased pseudospectral methods. All runs 
freely decay in time, with no external forcing. 

The type labelled RAN are 128 3 incompressible runs, 
with random broadband initial conditions. Four cases 
arc distinguished by their values of a m and o~ c , spanning 
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TABLE I: Parameters in the MHD simulations shown in the 
figures. RAN, OT, and ABC are described in the text. N 
is the resolution, v and 77 are respectively the kinematic vis- 
cosity and magnetic diffusivity, and a c and a m reffer to the 
normalized cross and magnetic helicities defined in the text. 
For the 2D run, o m is based on the mean-square flux function. 



Run 


N 3 


v — 77 




urn 


RANI 


128 3 


2.5 x KT 3 








RAN 2 


128 3 


2.5 x 10" 3 


0.5 





RAN3 


128 3 


2.5 x 1(T 3 





0.5 


RAN4 


128 3 


2.5 x 10" 3 


0.5 


0.5 


OT1 


128 3 


5 x 10" 3 


0.4 





OT2 


256 3 


1.5 x 1(T 3 


0.4 





OT3 


512 3 


7.5 x 10~ 4 


0.4 





ABC1 


128 3 


3 x 10" 3 





0.5 


ABC2 


256 3 


1.25 x 10" 3 





0.5 


ABC3 


512 3 


6 x 10~ 4 





0.5 


ABC4 


1536 3 


2 x 10" 4 





0.5 


2D 


1024 2 


2.5 x 10~ 4 









a range of possibilities for relaxation starting from a fully 
random state. 

We also employ two other types of initial conditions 
in which helicities are controlled and the fields are more 
ordered. OT runs are a generalization of the 2D Orszag- 
Tang (OT) vortex @, a standard large scale initial con- 
dition for MHD turbulence. In our OT case, initially 
energies E v = Ef, = 2, a c ps 0.4, and a m w 0. An- 
other set or runs labelled ABC consists of a parameter- 
ized large scale helical flow, an uncorrclated and helical 
large scale magnetic field, and added noise with energy 
spectra ~ fc~ 3 exp[-2(fc/fc )] 2 at t = 0, with k = N/6 
jig ]. These runs have E v = E^ = 0.5, a c ps 1 x 10~ 4 , and 
a m « 0.5, while numerical resolution and Reynolds num- 
bers vary (see Table I) . Finally we analyze a small spatial 
region near a current sheet in a very high Reynolds num- 
ber ABC simulation, ABC4 in the table. 

Probability density functions. Our main diagnostics are 
probability density functions (pdfs) of the local cosine of 
the angle 9 between v and b 

v ■ b 

cos 9 = cos(v, b) = - — — — (3) 
l v ll b l 

which are computed for each run. 

The distribution function for RAN2 is shown FigfTJ), 
at times t = 0,0.5, 1.0 and 2.0. These distributions are 
highly peaked near cos 9 ps 1, much more so than would 
be needed to account for the cross helicity which is ini- 
tially a c ~ 0.5, decaying to a c = 0.24 at t = 2.0. The 
more peaked curves are for the progressively later times. 
The results for RANI, having no helicities, arc shown 
in Fig. [21 Now, the distributions are suppressed near 
cos9 w and strongly peaked near cos# ps ±1 indi- 
cating an enhanced probability of magnetic and velocity 
field being very nearly aligned or antialigned. Enhanced 
directional alignment occurs even when the globally av- 
eraged cross helicity is approximately zero. We do not 



show results for RAN3 and RAN4, with a m k, 0.5, as 
the distributions are almost indistinguishable from the 
corresponding case with er rn ps 0. 
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FIG. 1: Pdfs of cos# for initial normalized cross helicity a c = 
0.5 for Run RAN2. Global normalized cross helicity is 0.24 
at t=2. Different lines are for different times (see text). 
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FIG. 2: Pdfs of cos 9 at times t=0 (dotted), 0.5 (dash-dotted), 
1.0 (dashed), 2 (solid) in a 3D simulation a c ~ and a m ~ 
(Run RANI). The initial distribution is flat. 

The pdfs in the OT runs (not shown) are asymmet- 
ric and strongly peaked at cos 9 pa 1, as in the RAN2 
and RAN4 cases. For the ABC runs, with no net cross- 
hclicity, the pdfs peak at cos# ~ ±1 after less than half 
a turnover time, following the pattern of the RAN runs. 
Figure [3] shows the pdfs from the ABC runs at the peak 
of dissipation (t ps 4) for different Reynolds numbers. 

This local alignment process is fast in all cases, with 
substantial and apparently nearly saturated alignment 
occurring in less than one large scale turnover time. As 
stated above, no clear dependence with the Reynolds 
numbers is seen when we compare cases ABC1-4. 

When pdfs of cos 9 are computed in the vicinity of a 
cluster of strong current sheets, or in regions of strong 
shear in the magnetic field (run ABC4), an only slightly 
different result is obtained (Fig. 2]). Inside the current 
sheet, the magnetic and velocity field are strongly an- 
tialigned (which gives the peak near —1), and the pdf 
is linear. As larger subvolumcs surrounding the cur- 
rent sheet are considered, or at later times when current 
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FIG. 3: Pdfs of cos6< in runs ABC1 (solid), ABC2 (dotted), 
and ABC3 (dashed) at the peak of energy dissipation and 
ABC4 (dash-dotted). 
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FIG. 4: Pdfs of cos# in the vicinity of a current sheet (sub- 
volume of 150 grid points) in a 1536 3 simulation with ABC 
plus noise initial conditions (run ABC4). 

sheets accumulate and interact, and thus more current 
sheets with different alignments are integrated in the sub- 
volume, the pdf converges towards the form seen in Fig. 
Q for the whole flow. 

Physics of alignment. Why does local alignment take 
place in these simulations? And why is it so fast? Manip- 
ulating the MHD equations, Eqs. ([1]) and ([2]) in the ideal 
case [y = ?/ = 0), one finds the equation for evolution of 
the local cross helicity: 

d(v ■ b) v 2 

\ + v • V(v ■ b) = b • V b VV. (4) 

dt y ' 2 w 

The terms on the left are the convective derivative, indi- 
cating that v • b is advected by the velocity field. The 
terms on the right are divergences: using that V • b = 0, 
and when integrated over volume with the proper bound- 
ary conditions (e.g. periodic boundaries) they vanish. 
This expresses the simple fact that the total cross helic- 
ity is an ideal invariant in MHD. 

However, gradients of kinetic energy and pressure gra- 
dients affect the local alignment between the two fields. 
The first term on the right of Eq. ((4]) shows that gra- 
dients in the kinetic energy (e.g., shear) can change the 
alignment between the fields when they are parallel to the 
magnetic field lines. Indeed, a magnetic field line (which 
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FIG. 5: Pdf of cos(v, b) from 30 years of Omni data (ISEE, 
IMP and other satellite data). Also shown is the pdf of 
cos(v, b) from Ulysses spacecraft data between 50 and 59 ° 
heliospheric North latitude during a polar pass in solar mini- 
mum conditions. 



behaves as a material line as follows from Alfven's theo- 
rem) tends to be distorted by the shear, and aligned with 
the velocity field locally. For a planar shear, this would be 
very similar to what is called field-line stretching. Pres- 
sure gradients aligned along magnetic field lines, from 
the second term on the right, also change the alignment. 
Where a pressure gradient is present, velocity goes from 
the region of higher pressure to the region of lower pres- 
sure. If the pressure gradient has a projection onto the 
magnetic field, the resulting velocity field will be aligned 
as a result with the magnetic field. Moreover, from di- 
mensional analysis we can estimate the time for the local 
alignment to take place as ~ bi/l, where bi is the ampli- 
tude of the magnetic fluctuations at scale I. 

Note that the induction equation is formally equiva- 
lent to the vorticity equation in hydrodynamics. Conse- 
quently by the same reasoning, it can be shown that the 
hydrodynamic velocity and vorticity fields tend to align 
locally, as found numerically in 0] for regions of low dis- 
sipation. This, replacing b by the vorticity u> — V x v 
in Eq. , occurs according to alignment of oj with gra- 
dients of the kinetic energy and the pressure. 

Solar wind observations. Using samples of spacecraft 
data we computed distributions of the alignment angle 
for two interplanetary datasets - the Omni dataset at 
1AU near Earth orbit in the ecliptic plane, and a sample 
of Ulysses data from high heliographic latitude. Fig. © 
shows the results of these analyses. The low latitude 
OMNI analysis is divided into intervals in which the large 
scale interplanetary magnetic field is directed either away 
from or towards the sun. The net Alfvenicity is outward 
at the higher latitude of the Ulysses sample. In each of 
these cases, the pdfs of the local alignment are consistent 
with the net cross helicity in each sample. 

Discussion and conclusions. The characteristic pdfs 
of cos 6 described above cannot be explained as a su- 
perposition of two uncorrclatcd Gaussian distributions 
for the velocity and magnetic fields, although the pdfs 
of the velocity and magnetic field themselves are Gaus- 
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FIG. 6: Cross-helicity density in a 2D incompressible MHD 
simulation, showing areas that have values of cos(v, b) < —0.7 
(black), |cos(v,b)| < 0.7 (gray) and cos(v, b) > 0.7 (white). 
Areas having highly aligned or anti-aligned velocity and mag- 
netic field fluctuations dominate the picture. 

sian (but clearly correlated). Pdfs computed from ran- 
dom broadband uncorrclated Gaussian-component veloc- 
ity and magnetic fields have a flat /(cos#) distribution. 
For the coherent ABC Hows, cos 9 peaks at initially, 
while for the non-helical RANI and RAN3 flows the dis- 
tribution is initially flat. All cases evolve towards the 
characteristic shape that is high-peaked at |cos#| = 1. 
In contrast, prior studies have shown that the distribu- 
tion of the induced emf, v x b, is accurately computed 
from the Gaussian statistics, for both high and low cross 
helicity, in simulations and in solar wind data [Til ] . What 
apparently accounts for the difference is that the induced 
emf does not correspond to a conserved quantity, while 



the alignment angle is closely associated with the ideally 
conserved cross helicity. The emf can be accounted for 
using Gaussian statistics, but alignment, even of Gaus- 
sian fields, is a dynamical quantity constrained by the 
local transport and conservation, as implied by Eq. (4). 

Note that Alfven vortices [12], which are coherent 
structures predicted for reduced MHD, have been re- 
cently observed in space plasmas [T3|; the generalized 
Alfven condition obeyed by these vortices corresponds 
to a local directional alignment. Evidently this type of 
robust alignment process may be influential in a variety 
of space and astrophysical plasmas in which turbulent 
relaxation operates, as well as in the neutral fluid case. 

We conclude that directional alignment is a rapid and 
robust process in turbulence. The magnetic and velocity 
fields respond, as described above, according to the lo- 
cal values of the shear and pressure gradients, essentially 
independently of the conditions at remote locations, lead- 
ing to local alignment or ant i- alignment; Fig. [^illustrates 
this localization or patchiness of the directional align- 
ment, using a 2D MHD simulation (see also [3). Since 
the alignment appears to be an essentially universal and 
rapid process, it would not be surprising if the coherent 
small scale structures in turbulence are associated with 
it. Indeed, the case shown in Fig. 3] is such an exam- 
ple where current sheets are observed to have maximum 
alig nment between the velocity and the magnetic fields 
[101 ]; similarly the local v-o; alignment may explain the 
slow return to full isotropy in fluid turbulence. 
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